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luid-induced instability problems, such as oil whirl and 
ail whip, still often occur as an unpleasant surprise in 
the field, but machines are regularly purchased with no 
stability test. 
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Ey far, the best way to test for stability of a rotating machine 
is to run it at its regular operating speeds and conditions. 
While at these conditions, perturb the rotor in the forward 
direction (with rotative speed) by a forward circular force (or 
moment) of known amplitude and phase at frequencies from 
slow roll to well over rotative speed, Observe the amplitude 
and phase of the resultant rotor motion at each perturbation 
speed. 

At each speed, the known input force vector at the pertur- 
bation speed is divided by the observed motion vector of that 
speed in order to yield the dynamic stiffness of the machine 
while at its usual operating conditions. 

The general dynamic stiffness equations of a simple rotor 
system (ignoring gyroscopics and assuming symmetry) are: 

Kpm = Ko- (w - AQ)? *M,, -w*Mp, 

Kouan = jDp * (w- a2) + {Dew 

where Kore i$ the Direct Dynamic Stiffness 

Kouan is the Quadrature Dynamic Stiffness 
Ko is the direct system stiffness 

w is the perturbation (or whirling) speed 

{118 the rotative speed 

A is the ratio of circumferential fluid whirling to 
rotative speed 

Moy is the fluidic inertial effect (not treated in 
this article) 
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M, 1s the effective rotor mass 

Dp is the bearing (or seal} damping 

Ds is the system damping away from the bearing 

(or seal) 

W is the total rotor weight 

The Direct Dynamic Stiffness and the Quadrature 

Dynamic Stiffness graphs are plotted separately as a function 
of the perturbation speed. The typical results are shown in 
Figure 1. 
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Figure 1 
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Da + Dg 
more complicated, but fundamentally the same, for typical rotar systemts, 
See the papers of Donald E. Bently and Or Agnes Muszynska for details. 


* More correctly w = in the example and a [ttle 


These direct and quadrature stiffnesses of this very simple 
rotor apply quite well, no matter how complicated the actual 
rotor system. For any rotative speed 9, the speed where the 
Direct Dynamic Stiffness Kog is zero is the well-known 
lateral resonance “critical” speed, point Q} in Figure 1. 

For that same rotative speed ©, there is also a zero of the 
Quadrature Dynamic Stiffness terms Kouan This is the fluid 
resonance, w = Ag", point @)in Figure 1. This fluid resonance 
is never seen except by perturbation methods like those 
described above. 

When these two resonances occur at the same speed w, 
then the rotor system whirls in accordance with well-known 
stability laws which govern the behavior of rotor systems as 
well as for any other closed loop control system. 


6 ORBIT 


Since this full perturbation test has not been used in the 
specification for purchasing of machinery, here is a very 
simple test that can be done to give a vital portion of the 
results of the full test. This test is to use the weight of the rotor 
as a steady state (w = 0) perturbation force. 


To accomplish this test, the gap voltages of the shaft observ- 
ing AY probes at each end of the machine are observed while 
the machine is on the test stand in its normal position. Next, 
the machine under test is rotated to 45° left and the probe 
gaps are again observed. The machine is rotated 45° right and 
again the probe gaps are used. Of course, this has the effect of 
rotating the steady weight of the rotor in the bearings so that a 
gravity vector can be resolved. See Fisure 2. 
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Figure 2 
Call the measurement with rotation to the right A, and with 
rotation to the left B, and the weight of the rotor W. The 
“horizontal” force, F,, is Fy = W sin 45°. See Figure 3. 
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Figure 3 
The motion vector change Zy is Zy = Za — Zp 
The steady state (w=()) stiffness vector for any rotative 
speed & Is 
Koo Z Fu £08: 
Zy /Azu 


Fy is considered to be referenced at 0°. This Koao 15 
separated into two components: 
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Now, throw away the direct term as it will usually contain 
the fluidic inertial effect term, which acts as a negative spring 
at w = 0 and confuses the results. 

The remaining term is the steady state quadrature stiffness 
term. This term is Kouan wo = AD, the fluid support 
stiffness term. Plot this term on a Direct and Quadrature 
Stiffness graph, such as Figure 1. Add to this graph the line for 
rotative speed, &. Look up the self balance resonance, WRES, 
speed, and the “first critical” for this machine. Plot this on the 
zero stiffness line of the Direct Stiffness plot. 


Now, look up the effective weight of the rotor, Were (usu- 
ally about 60% of actual rotor weight, the result of distributed 
weight between bearings, as compared to a concentrated load 
between bearings). Calculate the system effective direct stiff- 
ness Kor = (wgesl Mp 


Plot this Kpg on the w=0 line of the Direct Dynamic 
stiffness Graph. You can now draw the negative parabola, 
Kor TT w M Epp on the Direct Dynamic stiffness Graph. 


As tor the Quadrature Stiffness Graph, it is always a 
straight line with a slope, D, for damping since the entire line 
i$ Kouap = (w-AQ)D. Rotative speed, &, is known, and you 
have already plotted Kauan at w= 0. 

This leaves the A ratio to be determined. With no anti- 
swirling, a fair bet is a little less than ‘5, say 0.48. (Beware, 
some pumps and some high pressure compressors have A 
ratios of anywhere from 0.5 to 1.2!) Plot this AD term on the 
Kouan—o line of the Quadrature Stiffness Graph and draw a 
straight line from the Kouan «=o fluid wedge support stiff- 
ness term through the Ag point. This is the Quadrature 
Stiffness Graph. 


If the manufacturer has installed antiswirl to enhance sta- 
bility, this data may be done with and without the antiswirl. 
This will exhibit and measure the corrective effect (reduction 
of the quadrature stiffness) of the stability control method. 


The bottom line is that, for a stable machine, the AD term 
on the Koyan graph had better be SMALLER than the wpes 
point on the Ko graph! When they are at the same fre- 
quency, instability begins. If the AD term is larger, then wpps 
moves to a matching larger term merely by increasing the 
orbital size of the whip due to the instability. (A bigger size of 
orbit means a bigger bearing and seal stiffness, therefore, the 
higher “critical” speed. 

It is obvious that this simple test takes some careful han- 
dling of alignment, so as not to introduce another steady term 
besides gravity, and requires the modal probes to reflect the 
motion from the probe locations to the centers of seals or 
bearings. It is also obvious that you can have the PREDIC- 
TION of stability or instability if that is important to you. W 


